Radical-ion-pair reactions are the biochemical equivalent of the optical double slit 

experiment 
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Radical-ion-pair reactions were recently shown to represent a rich biophysical laboratory for the 
application of quantum measurement theory methods and concepts. We here show that radical- 
ion-pair reactions essentially form a non-linear biochemical double slit interferometer. Quantum 
coherence effects are visible when "which-path" information is limited, and the incoherent limit is 
approached when measurement-induced decoherence sets in. Based on this analogy with the optical 
double slit experiment we derive and elaborate on the fundamental master equation of spin-selective 
radical-ion-pair reactions that covers the continuous range from complete incoherence to maximum 
singlet- triplet coherence. 



I. INTRODUCTION 

The double slit experiment is the archetypal system 
manifesting the core concepts of quantum physics. We 
here show that Nature has already designed biologically 
significant chemical reactions that act as a double slit 
interferometer. Spin-selective radical- ion-pair reactions 
are perhaps the only example in chemistry where spin 
degrees of freedom and their relatively small interaction 
energy can have a disproportionately large effect on the 
outcome of chemical reactions. Their study is at the core 
of spin chemistry [l| , by now a mature research field di- 
rectly related to photochemistry [2| and photosynthesis 
Radical- ion-pair reactions determine the late-stage 
dynamics in photosynthetic reaction centers [4|, [5[ , and 
are understood to underlie the avian compass mecha- 
nism, i.e. the biochemical compass used by migratory 
birds to navigate through the geomagnetic field [6l-[l2j. 
A deeper understanding of the quantum dynamics inher- 
ent in radical-ion-pair reactions will thus enhance or even 
fundamentally alter our understanding of these biological 
processes. Other biochemical systems studied from the 
quantum interference perspective include electron trans- 
fer pathways [13], and in a much broader context the 
coherent control of molecular dynamics Il3 an d the fem- 
tosecond control of chemical reactions jlMTlrjl]. 

Radical-ion pairs are biomolecular ions created by a 
charge transfer from a photo-excited D* A donor- acceptor 
molecular dyad DA, schematically described by the reac- 
tion DA -)> D*A -» D* + A*~, where the two dots rep- 
resent the two unpaired electrons. The magnetic nuclei 
of the donor and acceptor molecules couple to the two 
electrons via the hyperfine interaction, leading to singlet- 
triplet mixing, i.e. a coherent oscillation of the spin state 
of the electrons. Charge recombination terminates the re- 
action and leads to the formation of the neutral reaction 
products. It is angular momentum conservation at this 
step that empowers the molecule's spin degrees of free- 
dom to determine the reaction's fate: only singlet state 
radical-ion pairs can recombine to reform the neutral DA 
molecules, whereas triplet radical-ion pairs recombine to 
a different metastable triplet neutral product. 

Theoretically, the fate of radical-ion-pair reactions is 



accounted for by the time evolution of p, the density ma- 
trix describing the spin state of the molecule's two elec- 
trons and any number of magnetic nuclei. It was recently 
shown [17j that quantum measurement dynamics are cen- 
tral in these biochemical reactions. A new master equa- 
tion was derived [17] based on quantum measurement 
theory, as the radical-ion-pair recombination process was 
interpreted to be a continuous quantum measurement of 
the spin state of the pair's electrons. This master equa- 
tion accounts for the spin decoherence of unrecombined 
radical-ion pairs. The kinetics of the recombination pro- 
cess, i.e. the loss of radical-ion pairs due to the forma- 
tion of neutral products, must also be taken into account. 
The treatment of this problem in [17j applies only to the 
case of maximal singlet-triplet coherence. Describing the 
continuous range of partial coherence down to the other 
extreme of maximal incoherence is what we need in order 
to provide a complete description of this system. We will 
here develop exactly this complete theoretical descrip- 
tion of radical-ion-pair reactions, shown to be an almost 
perfect analogy of Young's double slit experiment with 
partial "which-path" information. In Section II we reca- 
pitulate the open-quantum-system approach to radical- 
ion-pair reactions, elucidating the physics of the unavoid- 
ably present decoherence due to the intra-molecule mea- 
surement dynamics introduced in [17j and subsequently 
studied also in [18|. In Section III we introduce the re- 
action terms in the two extremes of maximum coherence 
and maximum incoherence. We then study the general 
case of partial coherence that naturally leads to the full 
master equation presented in Section IV along with some 
of its predictions in very simple cases. Finally in Section 
V we elaborate in detail on several aspects of the analogy 
with the photon's double-slit experiment. 



II. RADICAL-ION PAIRS AS AN OPEN 
QUANTUM SYSTEM 

The reaction dynamics and energy levels of radical-ion 
pairs are depicted in FigJTJ We neglect effects like diffu- 
sion, collisions, spin relaxation etc and we only consider 
the fundamental quantum dynamics inherent in the re- 
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FIG. 1: (Color online) (a) Energy level structure of radical- 
ion-pair reaction dynamics. A photon excites the singlet neu- 
tral precursor molecule DA into D*A, and a charge trans- 
fer creates the radical-ion pair. The excited vibrational lev- 
els (DA)* of the neutral DA molecule form the measurement 
reservoir, which (i) acts as a measurement device for the radi- 
cal pair's spin state, and (ii) acts as a sink of radical-ion pairs, 
i.e. in the event of recombination, the electron tunnels into 
a reservoir state and a fast spontaneous decay results in the 
ground state DA (which is the singlet product) and a photon 
emission. Similar for the triplet reservoir. 



combination process. In all experimental studies we have 
a macroscopic number of radical-ion pairs. Each one of 
them is a single open quantum system regarding the spin 
degrees of freedom, because the recombination dynam- 
ics inherent in each molecule disturbs a what would be 
a unitary spin evolution. Moreover, due to recombina- 
tion, radical-ion pairs form neutral chemical products, 
i.e. we deal with an ensemble with a changing number 
of open quantum systems, with these two aspects being 
intertwined. By properly identifying the system and the 
reservoir degrees of freedom in radical- ion pairs, we re- 
cently derived [17] a trace-preserving master equation de- 
scribing the evolution of the radical-ion-pair's spin state 
until it recombines: 



dp n 



dt 



ks + k T 



{pQs + Qsp-2Q s pQs) (1) 



Here T~L is the Hamiltonian embodying the magnetic in- 
teractions within the molecule, Qs is the singlet projec- 
tion operator (the triplet projection operator is denoted 
by Qt), and ks and kr are the singlet and triplet re- 
combination rates, respectively. Equation (pQ) embod- 
ies the decoherence or dephasing of the singlet-triplet 
coherence brought about by the intra-molecule contin- 
uous measurement induced by the recombination pro- 
cess. Due to this internal decoherence process, coher- 
ent superpositions are gradually turned into incoherent 
mixtures. To explain this in more detail, we first note 
that the singlet and the triplet reservoir enter the dy- 
namics on an equal footing [17j , the corresponding Lind- 
blad terms being -ks(pQs + Qsp ~ ^QspQs)/^ and 
— kr^pQr + QTp — 2QtpQt)/2' Due to the completeness 
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FIG. 2: (Color online) (a) Single-molecule quantum trajec- 
tory illustrating quantum jumps interrupting a unitary state 
evolution driven by the magnetic Hamiltonian 7i. At some 
random time t\ the measurement outcome is Qs = 0, and the 
radical-pair's spin state is projected (red square) to the triplet 
state, from which the unitary evolution commences. At some 
later time £2 (black circle) the radical-ion pair charge recom- 
bines and disappears, (b) The recombination event changes 
the density matrix normalization from one to zero, (c) and (d) 
Feynman diagrams describing a singlet and a triplet recom- 
bination, (e) and (f) Feynman diagrams describing virtual 
tunneling to singlet or triplet reservoir states responsible for 
singlet- triplet decoherence. 



relation Qs + Qt — 1, the sum of these terms results in 
the second, dissipative term of (pQ). Hence both reservoirs 
essentially measure one and the same observable, Qs, the 
measurement result being Qs = or Qs = 1. The total 
measurement rate is (ks + &t)/2. This continuous mea- 
surement interrupts a what would be a unitary spin evo- 
lution driven by the magnetic Hamiltonian 1-L. In other 
words, this unitary evolution is interrupted by quantum 
jumps produced by the projectors Qs and Qt- One illus- 
trative example of such a single-molecule quantum tra- 
jectory is shown in Fig. 2a. At some random time the 
radical-pair's state jumps to the triplet state (Qs = 0), 
and at some later time the radical-pair recombines, hence 
from that point on (Q s ) = (Qt) = Tr{p} = (Fig.2b). 
The average of all such trajectories (without the recom- 
bination event) is embodied in the master equation (pp). 
The recombination event is described by Feynman dia- 
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grams shown in Figs. 2c and 2d. The radical-ion pair 
state is coupled (by tunneling) to the reservoir state, 
which then emits a photon (or phonon) and the neutral 
product (ground state of DA or (DA) T ). The dephas- 
ing produced by the quantum jumps is described by a 
diagram like the one in Figs. 2e and 2f. Here the radical- 
pair "momentarily" tunnels to a singlet (triplet) reservoir 
state, from which it tunnels back to a singlet (triplet) 
radical-ion-pair state and the evolution commences on 
from there. This coupling to and back from the reser- 
voir is taken into account by the second-order perturba- 
tion in the system-reservoir coupling that leads to (pQ), 
as shown in [17j. To recapitulate, the excited vibrational 
levels of the singlet DA molecule and the triplet (DA) T 
molecule form the singlet and triplet reservoirs, respec- 
tively. Virtual transitions to both singlet and triplet 
reservoir states are responsible for singlet-triplet (S-T) 
decoherence, whereas real transitions to singlet (triplet) 
reservoir states cause singlet (triplet) recombination. 

The state change of unrecombined radical- ion pairs due 
to the intra-molecule quantum measurement dynamics is 
one part of the whole picture. The other is the change 
of the ensemble density matrix due to radical-ion pairs 
leaving the ensemble, i.e. recombining (reacting) away. 
We now come to the description of the reaction terms. 



III. REACTION TERMS 

Suppose that at time t an ensemble of radical-ion pairs 
is described by p t . Sure enough, in the following time 
interval dt we will obtain dns = ksdtTr{pQs} singlet 
products and dnr = hrdtTrlpQ^ triplet products. In 
other words, at time t+dtwe will have dns+driT less rad- 
ical pairs. We emphasize that the measurement of these 
product molecules represents classical information, the 
acquisition of which or not cannot have any back-action 
on the state of unrecombined radical pairs, as much as the 
detection of photons on the observation screen beyond 
Young's double slit does not back-react on the quantum 
state of the rest of the photons flying through the dou- 
ble slit. The information on the product molecules is in 
principle conveyed by the detection of d ns photons (or 
phonons) as depicted in FigsJ2t and [2]i [19]. The ques- 
tion we will address is the following: how is p t consis- 
tently evolved into pt+dt — Pt + dpi In order to simplify 
and motivate the present discussion we consider just two 
states [20j, the singlet \S) and the triplet |T), i.e. we 
consider a two-dimensional Hilbert space. Our results 
are generally valid for any radical- ion-pair Hilbert space. 



A. Maximum incoherence extreme 

Suppose we know that at time t we have an incoherent 
singlet-triplet mixture, i.e. p t = \s\S)(S\ + \t\T)(T\, 
with Xs + At = 1, or equivalently, using the projec- 
tors Q s = \S)(S\ and Q T = \T)(T\ we can write 



Pt = ^sQs + AtQt- This is the maximum incoherence 
extreme, i.e. there is zero S-T coherence. We would 
then detect dns — ksXsdt singlet and dnr — krXrdt 
triplet products. Due to our knowledge of p t , we would 
also know for sure that the detected dns singlet and dnr 
triplet products must have originated from singlet state 
and triplet state radical- ion pairs, respectively. Thus we 
would write 

dpincoh = -ksdtQspQs ~ k T dtQ T pQ T (2) 

That is, we project out the singlet and independently 
the triplet part of p t by the reacted fraction of singlet 
(ksdt) and triplet (krdt) radical-ion pairs. Obviously 
Tr{dpi nC oh} = —dns — dnx, as it should be. 

B. Maximum coherence extreme 

Suppose, on the other hand, that at time t all radical- 
ion pairs are in the maximally coherent state = (\S)-\- 
\T))/y/2 [21]. Thus p t = |^| = (\S)(S\ + \T)(T\ + 
\S)(T\ + \T)(S\)/2. Since Tr{p t Q s } = Tr{p t Q T } = 1/2, 
we expect dns — ksdt/ 2 singlet and dnr = krdt/2 
triplet products. What information does p t convey about 
the possible precursors of the dns singlet and dnT triplet 
products? None. As the reacted radical-ion pairs cease 
to exist, so does the information about their particu- 
lar quantum state just prior to recombination. So now, 
in order to update pt we have to remove the complete 
single- molecule density matrix p t /Ti{pt} as many times 
as many products we measured, i.e. 

dp C oh = ~(dn s + dn T )p t / r Ti{p t } (3) 

This is a crucial point. The spin state of the singlet or the 
triplet neutral product will indeed be Q s ptQ s /Tr{p t Qs} 
and Qt PtQr /TY{p t Q T } , respectively, but that does not 
imply that the state of the precursor radical-ion pair was 
a pure singlet or a triplet. It was a coherent superpo- 
sition. In other words, coming back to the double slit 
analogy, observing a photon in the upper (lower) part 
of the observation screen does not imply that the photon 
had crossed the upper (lower) slit. That would obliterate 
the physical essence of quantum superpositions. Again, 
as in the previous case, Tr{dp co h} = —dns — dnx, as it 
should be. It is finally noted that the treatment of reac- 
tion kinetics in our first publication [17j applies to just 
this special case of maximum singlet-triplet coherence. 

C. General case: partial coherence 

We now have to address the general case of partial 
coherence that spans the region between the above two 
extremes. Suppose that p t is an arbitrary density ma- 
trix in the radical-pair's spin space. What we need is 
a measure of singlet-triplet coherence of p tl that will 
continuously span the intermediate region between zero 
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coherence and maximum coherence. Before deriving 
this measure, we briefly mention the analogous situa- 
tion in the optical double slit experiment. Photon in- 
terference with partially coherent light is now a textbook 
example [22], the intensity at the observation screen is 
J(r) = I 1 +I 2 + 2y/hhi(x 1 ,x 2 ), where ^ - (l^^)! 2 ) is 
the intensity at the observation point r that would result 
from only slit j being open, with j = 1,2, and 7(^1, #2) 
the complex first-order coherence function, with | | = 1 
for complete coherence, I7I = for complete incoherence 
and < I7I < 1 for partial coherence. 

The analogous measure of S-T coherence for radical- 
ion pairs is derived as follows. By multiplying p t from left 
and right by the unit operator, and replacing the latter 
by Qs+Qt = 1, we can write p t = Pss+Ptt+Pst+Pts, 
where p S s = QsPtQs, Ptt = QrptQr, Pst = QsptQr 
and pTS — QrptQs- The coherence of p t is obviously 
Pst + Pts- Since psr and prs are traceless (because 
QsQt = 0), we consider the trace of their product as 
the measure describing the "amount" of singlet-triplet 
coherence, i.e. the quantity Tr{psT/>rs}- We define the 
S-T coherence measure p co h as 



Pcoh 



TrjpsTPTs} 
Tr{pss}Tr{p TT } 



(4) 



For the maximal coherence extreme it is p co h = 1, while 
for completely incoherent mixtures we have p co h = 0, 
with all other values covering the intermediate partial 
coherence regime. These properties of p co h are proved in 
Appendix A. 

So now, the only statement we can make about the 
change in the density matrix dp consistent with the in- 
formation at hand is that 



dp r = (1 - PcohMPincoh + Pcoh^Pcoh 



(5) 



This is the reaction term for a general state p t having 
partial S-T coherence. In other words, for partial S-T 
coherence, we only know the probability that, for ex- 
ample, a singlet product came from a singlet precursor 
radical-pair, or from an S-T coherent radical-pair. The 
former is 1 — p co h and the latter is p co h- 



IV. MASTER EQUATION 

To arrive at the complete theory, we have to add 
the change of the density matrix due to the recombined 
radical-pairs, dp ri with the change of state of the unre- 
combined radical pairs, dp nr , leading to the master equa- 
tion 



k 7 



(pQs + QsP - 2Q s pQs) 



(1 - Pcoh)(k s Q s pQs + KtQtpQt) 
■ Pcoh (k s Ti{Q sP } + k T Ti{Q TP }) 



Tv{p} 



(6) 



This master equation embodies the fundamental quan- 
tum dynamics of spin-selective radical-ion-pair reactions. 
Since dTi{p nr } = 0, it is seen that the trace of p decays 
as it should, i.e. according to the number of recombined 
radical- ion pairs, dTr{p} = — dn$ — dnx- The singlet 
and triplet reaction yields are then Ys = / °° dns and 
Yt = J °°(inT, respectively. The traditional theory (2^ 
identically follows from (j6]) by setting p co h = for all 
times. That is, the traditional theory by design com- 
pletely suppresses the effect of radical-pair spin coher- 
ence. We shall show (Section V.D.) that this is done by 
changing the state of unrecombined radical-ion pairs in 
an unphysical way. 



A. Predictions of Master Equation 

The predictions of this master equation and their com- 
parison with the traditional paradigm is a rather long 
program to be undertaken elsewhere. We will here illus- 
trate just one quite striking example on the differences 
between the new and the previous theoretical under- 
standing (a more recent approach by Jones & Hore [l8[ 
offers similar predictions as the traditional theory). This 
example serves to unravel fundamental differences in the 
quantum mechanical understanding of radical-ion-pair 
reaction dynamics in a simple enough setting that will al- 
low the analogy with the optical double-slit experiment to 
clearly illuminate these differences. We thus consider the 
most basic case of no magnetic interactions and a single 
recombination channel, for example the singlet, i.e. we 
completely eliminate S-T mixing and we take kr = 0. We 
chose as initial state the coherent S-T superposition [2l[ 
l^i) = (l^) + l^))/v // 2 (a coherent superposition is chosen 
exactly because it is in this maximum coherence regime 
Pcoh = 1 that the new master equation significantly de- 
parts from the previous theory) . We plot the expectation 
values of the projectors Qs and Qt in Fig HI It is evident 
that the traditional master equation of spin chemistry 
predicts that half of the radical-ion pairs stay locked in 
the non-reacting triplet state, whereas the other half pro- 
duce singlet recombination products. It is as if at t = 
we make a global von-Neumann measurement of the en- 
tire population of radical-pairs, which are then projected 
into either the non-reacting triplet and stay there for- 
ever, or into the reacting singlet, with probability of 0.5 
for both cases. We will explain shortly how this behavior 
comes about. In contrast, the master equation (j6j) pre- 
dicts that 75% of the molecules will react (FigHJ). This 
is because what actually happens is a continuous weak 
measurement operating within individual radical-ion pair, 
resulting to spin decoherence. Before the latter sets in, 
radical-ion pairs react through the singlet channel, but as 
time progresses, the initial coherent state is transformed 
into an incoherent mixture, so the molecules are gradu- 
ally locked in the non-reacting triplet state. This process 
can be visualized with the diagram shown in FigHJ de- 
scribing single-molecule quantum trajectories. We will 
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FIG. 3: (Color online) Radical-ion-pair recombination dy- 
namics in the absence of singlet-triplet mixing, (a) Time 
evolution of Tr{Qsp} and (b) of Tr{Qxp} as predicted by 
© (solid line) and the traditional as well as the Jones-Hore 
theory (dashed line). 
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FIG. 4: (Color online) Single radical- ion-pair state evolution 
with no singlet-triplet mixing, a single recombination chan- 
nel (kr = 0) and a coherent initial state. A radical-ion pair 
starting out in the state po either recombines (probability p r ) 
or it does not (probability p nr = 1 — p r . If it does not re- 
combine it will either be projected to ps (probability qs), or 
projected to pr (probability qr) or else remain at po (proba- 
bility 1 — qs — <?t), and the same scenario unfolds from there 
on. 



trace the state evolution at time steps dt, as shown in 
FiglU Starting with the coherent state po = with 
= (IS) + \T))/y/2, the probability for singlet recom- 
bination is p r = ksdt(Qs) = ksdt/2. In a particular 
realization of this single molecule experiment this recom- 
bination takes place in the first time interval (0, dt), and 
the reaction terminates. In some other realization the 
radical-ion pair will not recombine in this first time in- 
terval dt. Non-recombination happens with probability 
p nr = 1 — p r . According to (pQ), spin coherence will be 



lost for non-recombining radical pairs. In single- molecule 
realizations, this happens as follows. Conditioned on 
the fact that the radical-ion pair does not recombine 
within the interval (0, dt), there are three possibilities 
for what can happen: (i) a singlet projection to the state 
p s = QsPoQs/Tr{poQs} takes place with probability 
Qs — ksdtTr{poQs}/2 = k$dt/4, (ii) a triplet projection 
to the state pr = Qt PoQt /TrjpoQx} takes place with 
probability qr = ksdtTr{poQT} /2 = ksdt/4, and (iii) 
the state remains po with probability qo = 1 — ksdt/2 
[23j]. The average of all three possibilities, for which 
Qs H~ Qt H~ Qo = 1, reproduces the master equation (pQ) 
of unrecombined radical pairs. If the radical-pair is pro- 
jected to the singlet state ps-, it will have (Qs) = 1, so it 
will recombine for sure at some point. If it is projected 
to the triplet state pt, it will remain there forever, as 
k T = 0. Finally, if the state remains po, the same sce- 
nario will unfold in the next time interval. Since those 
molecules that have been projected to ps will definitely 
recombine at some point, it is easy to compute the sin- 
glet yield Ys by summing the total probability for re- 
combining from po with the total probability for being 
projected to ps- The probability of the state remaining 
po is po — p nr qo = (1 — ksdt/2) 2 w 1 — ksdt, whereas 
the probability for singlet projection is ps = PnrQs — Qs- 
Therefore 

Ys = Pr+PS J rPo{pr + Ps) + P^Pr + Ps) + 
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(1 - k s dt) 2 -k s dt - 



(7) 



This is how the result shown in Fig|3] is produced, i.e. 
this is how the total singlet yield comes out to be larger 
than 0.5, and the fraction of radical- ion pairs remain- 
ing locked in the non-reacting triplet state is accordingly 
smaller than 0.5. At first glance this result sounds unex- 
pected at the least. How can this be acceptable, since 
we start with 1/2 probability of the radical-pair spin 
state being singlet or triplet? According to the tradi- 
tional theory [25], coherent superpositions of the form 
ai|D #+ A #- )+a 2 |DA), where |D #+ A #_ ) is the radical-ion 
pair state and |DA) the neutral product state, are permis- 
sible. So the initial state (\S) + \T))/y/2 will evolve to 
|singlet product) + (triplet radical pair))/v / 2, hence the 
expected singlet yield is 0.5, equal to the final fraction 
of unrecombined triplet molecules. In reality, however 
[26| , there cannot be any quantum superpositions be- 
tween react ants and reaction products of an irreversible 
exergonic reaction. In other words, there is no physical 
operator having the radical-ion-pair state and the neutral 
product state as its eigenvectors. Put differently, there is 
no physical mechanism performing a global von Neuman 
measurement that would result into a 50-50 distribution. 
Hence there is no contradiction with projective measure- 
ment quantum mechanics of having a singlet yield larger 
than 0.5, as is actually the case in this example. In still 
other words, the intramolecule quantum measurement of 
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Qs takes place in the radical-pair's spin space and leads 
to singlet or triplet projections through Qs or Qt- These 
quantum jumps are completely independent of the recom- 
bination events, which can take place at any time, so that 
the reaction yields cannot be in any physical way consid- 
ered to represent the result of any kind of measurement 



B. Energy and Angular Momentum Conservation 

It might appear that the state of affairs previously out- 
lined might lead to non-conservation of energy or angu- 
lar momentum. For example, consider adding a singlet- 
triplet energy splitting produced by a spin-exchange cou- 
pling of the form Jsi • S2. Then the energy difference 
of the singlet and triplet levels would equal J, the ini- 
tial energy of the coherent superposition would be J/2 
(defining as zero the energy of the singlet level) and it 
would appear that the final energy balance is broken, 
since the final triplet population is 25%, i.e. an energy 
of 0.25 J/2 per molecule is missing. The same argument 
would apply for angular momentum conservation, had we 
considered a superposition of the form (\S) + |T + ))/V% 
where T + = | tt) is the plus-one projection of the total 
electron spin. We will briefly show why there is actually 
no violation of these conservation rules. First of all, such 
superpositions cannot be produced in the first place by 
any physical magnetic Hamiltonian since the latter 
will always couple states satisfying energy and angular 
momentum conservation. Nevertheless, suppose that in 
some way such superpositions are produced. We intro- 
duce a more general definition of p C oh- 



PcohO) 



|((Tr{psT(t)pTs(t + r)}))| 
Tr{pss}Tr{p T T} 



(8) 



where ((..)) denotes a time average over r (with r be- 
ing larger than the inverse S-T energy difrerence and 
smaller than characteristic time-scale of the reaction) and 
Prs(t + t) = e~ lUr PTs(t)e znT . In the above mentioned 
case of a singlet-triplet splitting by J, the term inside 
the time-average becomes e~ lJt , which quickly rotates in 
the complex plane at an angular frequency J, making 
Pcoh ~ 0. In the p C oh = case we retrieve the traditional 
theory, which does conserve energy. Essentially, an S-T 
energy difference makes the two "parties" in the quan- 
tum superposition distinguishable, suppressing the effect 
of quantum interference. Similar considerations apply to 
angular momentum conservation. 



V. ANALOGY WITH THE DOUBLE-SLIT 
EXPERIMENT 

We will now further dwell on several points concern- 
ing the analogy with the photon's double slit experiment, 
schematically depicted in FigJSl This analogy will pro- 



vide a comprehensive understanding of the quantum dy- 
namics of spin-selective radical-ion-pair reactions. 



A. Creation of Singlet-Triplet Coherence 

In the optical double slit experiment coherent superpo- 
sitions are created by free space propagation (i.e. diffrac- 
tion) across the slits. In radical-ion pairs coherent S-T 
superpositions are created by the magnetic interactions 
embodied in the magnetic interaction Hamiltonian H re- 
sponsible for the unitary evolution of p in (j6]). This 
can be easily shown as follows. We already pointed out 
that the density matrix can be written as p = p + p, 
where p = QspQs + QtpQt is the incoherent and 
p = QspQt+QtpQs is the coherent part of p. The Lind- 
blad term Qsp + pQs — 2QspQs appearing in the mas- 
ter equation (pQ) is equal to p, i.e. as already explained, 
the measurement dynamics dissipate singlet-triplet co- 
herence. The unitary part of the master equation (pQ), 
— i[H, p] is the generator of S-T coherence. Indeed, since 
QspQt + QrpQs = p, it follows from (pQ) that 



dp 



-iQs[H,p]Q T -iQT[H,p]Qs- 



-p (9) 



The first term of the first commutator on the right- 
hand side of (|9]) is (apart from —i) QsW-pQt, which 
after inserting Qs + Qt — 1 between T~L and p, and 
taking into account that Q 2 S = Qs and Q\ = Qt, be- 
comes UsspST + UstPtt, where H a (3 = QaUQp with 
a, j3 = 5, T. All other terms in the commutators of (|9j) 
are treated similarly, leading to 



dp 
~dt 



P~ 



i(j-LTsPss - Pss^-Lst 

+ ^STpTT — pTTHTS 

+ UssPst — Pst'Htt 
+ HttPts — pTsHss^j 



(10) 



Evidently, even if initially p = 0, i.e. there is no S-T 
coherence, pss or ptt will be different from zero, since 
Tr{p S s} + Tr{p T T} = Tr{p} > , so S-T coherence is 
generated by %st and %ts- At the same time, it is 
dissipated at a rate (ks + &t)/2. 

For simplicity in our previous discussions, i.e. for elim- 
inating the unnecessary complications of including spin 
interactions we just assumed an initial coherent super- 
position state without specifying how it was created. 
This is perfectly adequate for the purposes of illustrat- 
ing the fundamental quantum dynamics of the system 
and the conceptual differences from the previous physi- 
cal paradigm of the traditional theory. 



B. What is it that interferes? 

In the general case of partial singlet-triplet coherence, 
a singlet or triplet reaction product can be traced back to 
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(b) 



|S>-^S 



radical-ion 
pairs# • • 
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42 — 



singlet 
products 



triplet 
products 



FIG. 5: (Color online) Analogy between (a) the photon dou- 
ble slit experiment and (b) spin-selective radical-ion-pair re- 
actions. The magnifying glass denotes a detector acquiring 
"which-path" information. In the radical- ion-pair case this 
detector is realized by the singlet and triplet reservoir states, 
which continuously "measure" the pair's spin state at a rate 
ks/2 and hrj% respectively. 



Which-path information and destruction of 
interference 



How is this interference suppressed? In the double slit, 
if there is a photon detector providing definite "which- 
path" information, the interference pattern will be de- 
stroyed. Incomplete "which-path" information will re- 
duce the fringe visibility. Whether a photon detector is 
placed at the upper or the lower slit is indifferent, what 
matters is the total measurement strength of the pho- 
ton's path. In other words, it is immaterial whether a 
strong measurement is done at the upper slit and a weak 
measurement at the lower slit or vice versa, it is the total 
measurement strength that provides "which-path" infor- 
mation. This is so because in the extreme of zero mea- 
surement strength at one slit, i.e. when using just one 
photon detector e.g. at the lower slit, the absence of a 
detection implies that the photon went through the up- 
per slit, so "which-path" information is still available. In 
radical-ion pairs, it is the measurement rates ks and kr 
that gradually provide " which-path" information and re- 
duce Pcoh- Along the reaction, the internal decoherence 
process (a single- molecule process), taking place at the 
rate (ks + hr)/2 will gradually set in and damp singlet- 
triplet coherence that is generated by the magnetic in- 
teractions. In this process it does not matter whether 
ks is small and kr large or vice versa. What matters is 
the total measurement rate (ks -\-kx)/2 of the observable 
Qs, and this is what enters in Eq. (pQ). 



two different origins: (i) a singlet or triplet radical pair 
and (ii) a singlet-triplet coherent radical pair. By detect- 
ing a singlet or triplet reaction product, it is physically 
impossible to differentiate (i) from (ii), and this is embod- 
ied in the update rule dp = (1 - p co h) dp i nC oh+ Pcohdp co h • 
This update rule can also be read as dp = <ipi nC oh + 
Pcoh(^Pcoh — ^Pincoh), making evident the analogy with 
the double slit, i.e. the density matrix changes by the 
first, incoherent term plus the second term which is pro- 
portional to the S-T coherence p co h- The change dp will 
update the density matrix p t , which will govern the cre- 
ation of the reaction products in the following time in- 
terval dt. 



However, a point where the analogy with the optical 
double slit is not one-to-one is the fact that in the radical- 
ion-pair case there is no equivalent of the " photon- field" , 
hence there is no spatially dependent interference pat- 
tern. The distribution of photon clicks on the observa- 
tion screen has its equivalent on just two numbers: the 
time-integrated reaction yields. They are given by the 
combined contribution of all possible ways to reach the 
final singlet or triplet product, as depicted in the example 
of FiggJ 



D. Interference is a one-particle effect 

Radical-ion pair recombination is a single molecule 
process, similar to the double slit quantum interference 
which is a single photon effect. In that case, the de- 
tection of a photon on the observation screen after the 
two slits does in no way affect the quantum state of 
the rest of the photons. However, such a statement 
is embedded in the traditional theory. Here it is why. 
Suppose that at time t we have N radical-ion pairs all 
prepared in the coherent state = (\S) + \T))/y/2. 
Hence at this very instant in time p = iVpi, where 
pi = (\S)(S\ + \T)(T\ + \S)(T\ + \T)(S\)/2 is the single- 
molecule density matrix. No matter what the mas- 
ter equation satisfied by p is, we can obviously write 
that the change of p in the following time interval dt 
is dp = pidN + Ndp\. The first term represents the 
change in the density matrix due to reaction and the sec- 
ond the state change of the unrecombined radical pairs. 
Let us compare the predictions of our approach against 
the traditional theory. Since (Qs)t = 1/2, both theories 
predict dns = Nksdt/2 singlet products during dt. We 
have the machinery to calculate dpi. Indeed, from (pQ) we 
find dpi = —(ksdt/2)pi^ co \ l1 and hence 

dp = -dn s pi - -^^Pcoh (11) 
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where p coh = Np hcoh and p ljCoh = (\S)(T\ + |T)(5|)/2 
is the coherence of p\. Equivalently, we can arrive at 
(pTj) starting from the general master equation ((6]) and 
noting that at this instant p co h = 1- The interpretation 
of (jTTJ) is that we lose dns radical- ion pairs and the rest 
(unrecombined) lose some of their coherence. In contrast, 
the traditional master equation predicts that 

dp=-dn sP i+N^(\T)(T\-\S){S\) (12) 

Here the second term is highly problematic. Identifying 
this term with Ndpi, it follows that the unrecombined 
molecules have become more triplet and less singlet. The 
analogous statement in the double slit experiment would 
be that acquiring "which-path" information with a pho- 
ton detector placed at the upper slit (only k$ ^ 0) and 
measuring a particular number of " clicks" on the observa- 
tion screen, renders the rest of the photons flying through 
the screen more "down" and less "up". This is impossi- 
ble. In reality, if "which-path" information is acquired, 
the photons originally in the state (|up) + |down))/V2 
loose their coherence, i.e. their state approaches the 
mixed state ^|up)(up| + ^ | down) (down |, and that's why 
the fringe visibility is lost. The equivalent situation in 
radical-pairs is the loss of spin coherence described in 

(ED- 

E. On the non-linearity of the master equation 

Finally, we comment on the non-linear nature of (|6]), 
which should not come as a surprise (HI. We saw 
that if the initial state is = (\S) + |T))/V% 25% 
of the molecules are finally locked in the non-reacting 
triplet state. For a second ensemble initially in the 
state ^2) = (\S) — \T))/y/2, we arrive at the same re- 
sult. Now imagine mixing these two ensembles. If we 
naively take for the resulting system p = p\ + /?2, we 
find p = + |T)(T|, i.e. the final ensemble appears 

to be an incoherent singlet-triplet mixture. Since only 
the singlet reacts, at the end we would expect to have 
p t=OQ = |T)(T|, i.e. 50% of the total number of radical- 
ion pairs do not react. So what is the final triplet pop- 
ulation, 25% or 50% ? Let us answer by asking a differ- 
ent question: did the fact that we added the contents of 
two boxes together destroy the singlet-triplet quantum 
coherence? Not really. The ensemble p pr0 p = pi + P2, 
specifically prepared as described previously, is a proper 
mixture [29-31], i.e. each and every molecule is in some 
coherent S-T superposition, either or \^2)- We just 
don't know which. But the reaction does. The reaction 
is a single-molecule process, not some global process act- 
ing on all radical-ion pairs altogether. So at the end, we 
will indeed be left with 25% of the molecules we started 
with, whether we mix the two ensembles or not. On the 
other hand, if we start with the two ensembles and 
^2) and apply a decoherence process, so that the coher- 
ent superpositions are turned into incoherent mixtures, 



we would again obtain pimprop = 1*5) (*5 1 + l^)(^l- This 
time however, we have an improper mixture, not due to 
state ignorance, but due to a genuine loss of coherence, 
and we would get 50% singlet products and 50% triplet 
radical pairs. In general, one must clearly define how the 
initial state of a radical-ion-pair ensemble is prepared in 
order not to arrive at inconsistencies between theoretical 
predictions and observations. 

Returning to the double slit analogy, the optical inter- 
ferometer is linear since the interference pattern is de- 
pendent on the phase of the light field. Indeed, suppose 
we do the experiment in two consecutive runs with equal 
number of photons in each run. If for the second run 
we insert a tt phase shifter after one of the slits, the in- 
terference pattern will shift, and adding both runs there 
will be no interference pattern. That does not mean that 
the photon coherence was destroyed. This is the ana- 
log of the proper mixture of radical-ion pairs, where the 
singlet- triplet coherence is actually there. However, in 
this case, the reaction is a non-linear interferometer, so 
a phase change of the initial state does not alter the re- 
action yields (dns, dnr and p co h = 1 are the same for 
both and 1^2) )• If on the other hand, we destroy 
the coherence (|up) + |down))/\/2 of the photon state, 
e.g. with a detector providing "which-path" informa- 
tion, the interference pattern will genuinely disappear. 
This is the analog of the improper mixture of radical-ion 
pairs, where some decoherence process (for example the 
unavoidable one internal in the molecules) irreversibly 
destroys the coherence. 

VI. CONCLUSIONS 

Concluding, we have elucidated the fundamental role 
spin coherence plays in determining the fate of radical- 
ion-pair reactions, which were shown to be a non-linear 
spin-coherence interferometer, analysed in direct analogy 
with the optical double slit experiment. We have derived 
the fundamental master equation describing these biolog- 
ically significant biochemical reactions, shown to be the 
analog of the double slit experiment with partial "which- 
path" information, stemming from the singlet-triplet de- 
coherence brought about by the intra-molecule quantum 
measurement dynamics. As long as there are additional 
spin relaxation mechanisms, that suppress the singlet- 
triplet coherence p C oh, or in general, as long as p co h is 
small, the traditional theory of spin chemistry is expected 
to be successful. Significant deviations are expected to 
occur in situations where p co h is appreciable. 

Appendix A: Properties of p co h 

For the maximum incoherence extreme we have p — 
QspQs + QtpQt, i.e. by definition p S r = Pts = 0, 
hence p co h = 0. A maximally coherent state can be writ- 
ten as p = with 1*0) = a\S) + /3\T), and it is easily 
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seen that TrjpsTPTs} = |a| 2 |/3| 2 , whereas Tr{p S s} = M 2 
and Tr{pTT} = \fi\ 2 ! , hence indeed p co h = 1. The gen- 
eral case is treated as follows. Since psr and prs are 
traceless, Tr{p} = Tr{p S s} + Tr{p T T}, thus Tr{p} 2 = 
Tr{p SS } 2 + Tr{p TT } 2 + 2Tr{pss}Tr{p TT }. Second, since 
the projectors Qs and Qt are orthogonal, it is readily 
shown that Tr{p 2 } = Tr{p| s } + Tr{fe T } + 2Tr{p S T/>rs}. 
But Tr{p 2 } < Tr{p} 2 , hence 



Tr {Pssl + Tr{p 2 T } + 2Tr{p S TPTs} 
Tr{p SS } 2 + Tr{p TT } 2 + 2Tr{pss}Tr{p TT } 



<1 (Al) 



The first two terms in the nominator are less or equal the 
corresponding first two terms of the denominator, hence 
the same must hold true for the third terms (otherwise 
the ratio would become larger than one in the case of a 
pure state where Tr{p| s } = Tr{pss} 2 and Trjp 2 ^} = 



Tr{pxT} 2 )- We have thus shown that 



< Pcoh 



TrjpsTPTs} 
Tr{p S s}Tr{p TT } 



< 1 



(A2) 



To motivate the definition of p C oh, we consider the case of 
mixing N molecules in the singlet state, described by p = 
| S) (S | with another N molecules in the coherent superpo- 
sition p = |^}<V>| = (\S)(S\ + \T)(T\ + \S)(T\ + \T)(S\)/2. 
Take ksdt = 1/2, so the probability for the former to re- 
combine within dt is 1/2, while the respective probability 
for the latter is half as much, i.e. 1/4. If we detect a sin- 
glet recombination product, it will have thus originated 
with probability 2/3 from the singlet molecules and with 
probability 1/3 from the S-T coherent molecules. The 
latter probability is indeed derived from the expression 
for p coh . 
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